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Abstract. We use the semi-classical approach to study the non-equilibrium dynamics of the 
0(3) non-linear sigma model. For a class of quenches defined in the text, we obtain the order- 
parameter dynamical correlator in the thermodynamic limit. In particular we predict quench- 
dependent relaxation times and correlation lengths. The approach developed for the 0(3) non- 
linear sigma model can also be applied to the transverse field Ising chain, where the semi-classical 
results can be directly compared to both the exact and the numerical ones, revealing the limits 
of the method. 
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1. Introduction 

The out-of-equilibrium physics of many-body systems has attracted a lot of interest in recent years, 
not least because of the experimental realizations of ultracold atomic gases in optical lattices 
[1, 2, 3, 4, 5, 6]. These experiments have observed the dynamics of many-body systems on a 
long time scale after a quantum quench, finding essentially unitary time-evolution. In three- 
dimensional systems fast relaxation towards a thermal steady state has been observed. On the 
contrary, in quasi-one-dimensional systems the relaxation process is normally much slower and 
leads to a peculiar non-thermal stationary state [2]. These results have led to a huge theoretical 
push [3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25] to address 
fundamental questions such as whether there is an asymptotic stationary state, and, if it exists, 
which ensemble characterizes it. The belief is that observables of non-integrable systems effectively 
thcrmalize, which implies that their stationary state is characterized by a thermal Gibbs ensemble. 
Numerical works on non-integrable systems confirm this expectation, even if some contradictory 
results point out that some issues have not been completely understood yet [19, 20, 21, 26]. On 
the other hand, in integrable systems, because of the existence of local integrals of motion, the 
stationary state is expected to be described by a generalized Gibbs ensemble (GGE), where each 
mode associated with a conseved quantity is characterized by its own temperature. So far results 
on integrable systems have been focusscd on free fcrmion models, such as the transverse field Ising 
chain and the quantum XY chain [27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44]. 
For the transverse field Ising chain in thermal equilibrium Sachdev and Young first introduced 
a semi-classical description of the physical properties of the model in terms of ballistically 
moving quasi-particles [45]. This approach turned out to be incredibly accurate in predicting the 
temperature dependence of correlation length, relaxation time and in general to compute the order- 
parameter two point function in the ferromagnetic and paramagnetic phases. For global quenches 
this technique has been used to describe the dynamics of the transverse field Ising chain and the 
quantum XY model, with great accuracy [30, 37, 39, 44]. In these works quantitative features of the 
relaxation process have been explained with a quasi-particle picture, which had been introduced 
before in particular to study the evolution of the entanglement entropy [32, 46, 47, 48]. In practice 
the quench injects an extensive amount of energy into the system, which creates quasi-particles 
homogeneously in space, that then move ballistically with constant velocity. Because of momentum 
conservation, these quasi-particles are created in pairs with opposite momenta and are quantum 
entangled (within each pair). Their dynamics can be treated classically as long as they do not 
collide. But since collisions are unavoidable in one dimension, every scattering process has to be 
treated quantum mechanically. We shall use this semi-classical approach to calculate the dynamical 
correlation functions analytically for the 0(3) non-linear sigma model. The dynamics of this model 
have already been studied at finite temperature in equilibrium [49, 50], where the order-parameter 
correlation function shows a universal form. Here we shall study its behaviour after a quantum 
quench, preparing the system in a state which is not an eigenstate of the Hamiltonian. In addition 
a section will be dedicated to the case of the transverse field Ising chain, whose order-parameter 
two-point correlation function can be derived straightforwardly from the general one obtained for 
the 0(3) non-linear sigma model. 

The paper is organized as follows: first we introduce the model and the general form of the initial 
states we shall consider throughout the whole article, before introducing and commenting the main 
result of the paper, concerning the 0(3) non-linear sigma model. Then we describe the basic ideas 
behind the semi-classical technique and see in detail how to apply them in the present case. 
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2. The Model 

Let us start by considering the one-dimensional 0(3) quantum rotor chain, which is described by 
the following Hamiltonian: 

H^or = Jlj2^-J^n i n i+1 , (1) 

i i 

where hi is the position operator of the rotor on site i = 1, . . . , L with the constraint hf = 1 Vj, and 
Li = hi x f>i is the angular momentum operator. J is an overall energy scale and g is a positive 
coupling constant. The operators which appear in the Hamiltonian satisfy the usual commutation 
relations [Lf 7 L^] ~ ie a p^L] and [Lf 7 h^] — ie a ^nj , where a,(3,j represent the three spatial 
directions. The continuum limit of this lattice model is the 0(3) non- linear sigma model (nlcrm), 
whose Lagrangian density reads: 

C=^{d,hi)\ h} = \, (2) 

where hi = hi{x, t) are three scalar fields and g is a (bare) coupling constant. Here we have already 
set the maximal propagation velocity c = 1. This model is 0(3)-symmetric, renormalizable and 
asymptotically free, and it has three massive particles in the 0(3)-multiplet. The exact iS-matrix of 
this model is known [51], with any scattering event involving no particle production and the general 
n-particle 5-matrix factorizes into a product of two-particle amplitudes. It is worth noticing that 
the 0(3) non-linear sigma model also provides a description of the low-energy excitations of the one- 
dimensional antiferromagnetic 5=1 Heisenberg chain [52] , whose dynamical correlation functions' 
lineshape can be measured experimentally. In the g 1 limit, in the ground state of Hamiltonian 
(1) all rotors must be in the L\ = state to minimize the kinetic energy. The low-energy excitations 
above the ground state form a triplet with quantum numbers L\ = (—1, 0, 1). It is worth remarking 
that the structure of the low-energy spectrum is the same for any g > and the system has a 
gap A(<?). A finite gap is a necessary condition to apply a semiclassical approximation, as we shall 
see in the next section (for a deeper introduction to the semi-classical method and its range of 
applicability see [50] and references therein). 

In this article we study the dynamical correlator of the order parameter h z after having prepared 
the system in a squeezed coherent state, namely!: 




where |0) is the ground state of the model, K ab (k) is the amplitude relative to the creation of a pair 
of particles with equal and opposite momenta, Z\(k) are creation operators of an excitation with 
quantum number a = (—1,0, 1) (the z-component of the angular momentum) and momentum k. 
The main reason for choosing such an initial state comes from its relation with boundary intcgrable 
states: as shown by [10, 47] some dynamical problems can be mapped into an equilibrium boundary 
problem defined in a strip geometry, where the initial state \ip) acts as a boundary condition. In 
integrable field theories the most natural boundary states preserve the integrability of the bulk 

X Here for a matter of simplicity we ignore the existence of zero-rapidity terms in the definition of \ip). In addition 
it is known that an intcgrable boundary state as written in equation (3) is typically not normalizable, because the 
amplitude K does not go to zero for large momenta. Therefore one has to introduce an extrapolation time in order 
to make the norm of the state finite, as was done for instance in [18]. 
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model, or in other words, do not spoil the integrals of motion. These boundary states were originally 
studied by Ghoshal and Zamolodchikov [53], and are supposed to capture the universal behaviour 
of all quantum quenches in integrable models. Furthermore in [18] it was shown that any quantum 
quench of an integrable field theory with this kind of initial states leads to a steady state which is 
described by a GGE ensemble. Precisely they show that long time limit of the one point function of 
a local operator can be described by a generalized Gibbs ensemble§. In short, their result is strong 
evidence that Rigol et al.'s conjecture does hold for integrable field theories. 

The amplitude K ab (k) in expression (3) is a regular function which must satisfy a set of constraints 
that depend on the 5-matrix, such as crossing equations, boundary unitarity and boundary Yang 
Baxter equation. Different solutions of these equations form the set of integrable boundary 
conditions of the theory. || 

The integrable boundary states belong to the class (3), but the semi-classical approach applies to 
the larger class of quenches whose initial states are expressed by a coherent superposition of particle 
pairs. The order parameter n z (x) may be written as [50]: 

n z (x,t) cx (Zl(x,t) + Z {x, + (4) 

where Zo(x,t) is the (time-dependent) Fourier transform of Zo(k) and the ellipses represent 
multiparticle creation or annihilation terms, which will be considered negligible because they are 
subdominat in the long-time limit. Relation (4) comes from the observation that the operator n z (x) 
either creates a quasiparticle (Lj = 1) with Lf = at position x with some velocity v or destroys 
one already present in the system, as will be discussed in the section (4). 

The quench protocol is the following: At time t = we prepare the system in the initial state (3) 
and for t > this state evolves according to the Hamiltonian (1): 

Mt))=exp(-i#t)M (5) 
On the other side the time-evolution of an operator is written as: 

6{x, t) = cxp (iHt)O(x) cxp (-iHt). (6) 
We are going to analyze the two-point correlator: 

C" 1 ™(x 1 ,t 1 ;x 2 ,t 2 ) = {Mn z {x 2 ,t 2 W{x 1 ,t i m- (7) 

The autocorrelation function is obtained for x\ = x 2 = x', whereas the equal-time correlation 
function is defined by imposing t\ = t 2 = t', even though in this paper the emphasis will be given 
to the non-equal time case. 



3. Summary and discussion of the results 

In this section we summarize the main result of the paper. The semi-classical method is based upon 
the existence of a small parameter, namely the average density of excitation pairs with quantum 
numbers (a,b) n ab (k), in the initial state \if>): 

<«« SM -W, (8) 

§ A rigorous proof that the LcClair-Mussardo series work for the one-point functions in the case of a quench has 
been obtained in [54]. 

|| For instance for the case of the 0(3) non-linear sigma model with free and fixed boundary conditions these states 
have been studied in [55]. 
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which is valid in the limit max*; \K ab {k)\ 2 <C 1, Va, b. This limit also defines what we call small 
quench. The explicit computation of (8) will be carried out in appendix A. In the limit of small 
quenches we obtain the following expressions for the (non-equal time) two-point function of the 
order-parameter of the 0(3) non-linear sigma model (7) in the limit T = t 2 + t\ — > oo: 

C nlCTm (ar; t) = C propag (x, t) R(x; t), (9) 

where Cp Iopag (x 7 t) corresponds to the coherent propagation of a quasipaticle, while R(x;t) is the 
relaxation function, which describes the scattering with the excited quasiparticlcs. By definition 
t = t<i — t\ and x = X2 — X\ (without any loss of generality we assume both (x,t) to be positive). 
The pre-factor in (9) reads: 

C propag (a;,t) cx K {A^/x 2 - t 2 ), (10) 

where Kq is the modified Bessel function of the second kind and A is the mass-gap. The relaxation 
function R(x; t) is given by: 

/TT (^(k / POO \ 

— exp i-t(l - cos (j>)J — f k v k O[v k t- x]J 

/ [°° dk \ ( f°° dk \ 

x exp ( — x(l — cos0) J — fk Q[x — Vkt]j cos I xsincj) J — fk) (11) 

1 + 2cos(^)(P°° - 1)Q + cos(2(/>)(Q 2 - P°°) - P 00 Q 2 
X 1 - 2P 00 cos(20) + (P 00 ) 2 

Q is defined by Q — X)a=±i o ^° A ' /* = b fk b i v k is the velocity of the particle with momentum 
k and P ab = f °° |f P ab (k), where 

P ab (k) = - /oo fc „ ■ (12) 

El aK f ab 

Jo 

a.b 

Formula (11) is the main result of this article, and it describes the long-time behaviour of any initial 

state of the form (3). Nonetheless if we want the state to respect the 0(3)-symmetry we must 

impose P 00 = P 1 ' -1 = p^ 1 ' 1 = 1/3 (Q = 1/3), and all other probabilities equal to zero. 

The most general expression for the relaxation function R(x; t, T) for generic values of the times 

(ti)t2) will be given in the appendix B. Formula (11) defines the quench-specific time and length 

scales: 

i f°° dk „ „ , f°° dk , 

r- 1 = / —v k f k r 1 = / —fk- (13) 

Jo T Jo K 

We also note that the correlator of equation (11) is never thermal unless one prepares the system 
at t = in a thermal inital state. This result confirms the belief that one-dimensional integrable 
systems relax towards a peculiar non-thermal distribution, namely the generalized Gibbs ensemble 
(GGE). 
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4. The Semi-classical Theory 

The main idea behind the semi-classical approach to the computation of out-of-the-equilibrium 
correlators is basically encoded in the following representation of the correlator (7): 



where the function P({x v , k v , is the probability density of having quasiparticles at postion x v 
at time t — with momentum k v and quantum number A„. The matrix element in the equation 
represents the value of the correlator once we specify the particular configuration {x v , k v ,X v }. N in 
(7) is the total number of particles, and of course in the above equation an average over TV should 
be taken. The summation over {A„} averages over all possible quantum numbers of the N particles, 
which are the three possible values of Lf = (— 1, 0, 1). We will start by assuming the system to be 
finite with free boundary conditions, and only later will we take the thermodynamic limit. In spite 
of this we will assume the system to be translationally invariant, corrispondingly corrections due 
to the presence of boundaries will be considered negligible (for large enough systems). Therefore 
the correlator of equation (7) becomes a function of \x2 — x\\, whilst the same cannot be said for 
the time dependence, as was true in the equilibrium case. 

In equation (14) we have substituted a complicated time-dependent matrix element with a sum over 
all possibile initial states, which can be represented by pairs of ballistically-moving quasiparticles. 
This technique has already been used extensively to compute finite-temperature correlators (see for 
example [49, 50, 56]), and dynamical correlation functions [37, 39]. The general idea is the following: 
each quasiparticle carries a momentum k and a quantum number A„, and it is created at time t = 
together with a partner with equal and opposite momentum. If the squared modulus of Fourier 
transform \K ab (x)\ 2 of the amplitude K ab {k) is a fast enough decreasing function as |x| — » oo, the 
probability of creating these quasiparticles far apart from each other becomes negligible. For this 
reason in the follow we will always assume quasiparticles to be created in pairs in the same position 
x v , where the index v labels the quasiparticle pairs. 

We start by considering this gas of quasiparticles to be very dilute by tuning the amplitudes K ab (k), 
therefore the quasiparticle are entangled only within each single Cooper pair. This means that the 
probability of a particular initial state {x v ,k v ,X v } can be factorized in a product of single-pair 
probabilities: 



where in the last equality we made use that the system is homogeneous in space (corrections will be 
present close to the boundaries, but they are negligible in the thermodynmic limit). The quantity 
j^P(k v , X v ) is the probability for a single pair to be created at a certain position x v with momenta 
(k u ,~k u ) and quantum numbers (X V ,X' V ), which is graphically represented in figure (1). Without 
any loss of generality we can suppose x\ = 0, thanks to the (quasi)-translational invariance of the 
system. 

Calculating the matrix elements for the operator h z (oc cos 9) with the first few sherical harmonics, 




(14) 



M 



P({x v ,k v ,X v }) = \\_P{x v ,k v ,X v ) 



1 

- JmX\ P ^X u ) 



(15) 
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Figure 1. Graphical representation of the probability function P ab (k). 



one finds that either creates a quasiparticle with L\ = A = (L, = 1) at x = x\ and t = t\ with 
some velocity v or destroys one already present in the initial state {x u , k v , X u }. We assume the latter 
to be negligible, because the excitations are very dilute in space and therefore their density is very 
small (in other words we will use K ab (k) as an expansion parameter). This adjoint particle can be 
created by h z (xi,ti) either inbetween two different quasiparticle pairs or within two quasiparticle 
that belong to the same Cooper pair (i.e., two quasiparticles that originated in the same point x v at 
time t — 0). A possible dynamical scenario is pictured in figure (2). Due to the collisions with the 
other excited particles there arc only certain configurations where the quantum mechanical overlap 
in (tjj\h z (x2, t 2 )n z (xi, ti)\4>) will be non-zero. Similar to other models the 0(3) non-linear sigma 
model in the long-wavelength limit has a purely reflective scattering matrix: 

SaIaI (-l)^K 1 - (16) 
This means that in this limit the particles arc impenetrable and the sequence of {A„} does not change 
in time (see figure (1)). Here we are implicitly assuming that the momentum distribution fj. is a 
peaked function around zero. The operator h z (xi,ti) creates a quasiparticle with a probability 
amplitude eo(w), where the subscript o refers to the quantum number of the created particle and 
v is its velocity. This particle, together with the other quasiparticles, propagates under the action 
of exp (—iHt) and collides with them (in one-dimensional systems collisions amongst particles can 
never be ignored). The 5-matrix takes on an exchange form (16), and therefore particles only 
exchange their velocities while conserving their internal quantum numbers (see figure (2)). 
As a consequence, at any time t > t\ precisely one of the particles will have the velocity v of the 
particle which was created by n z (xi,ii), and will be at position x\ + v(t — t\). This very particle 
must be annihilated at time i 2 by h z (x 2 , t 2 ), otherwise the final state after the action of exp (iHt) 
will be orthogonal to the initial one. The probability amplitude that this particle is annihilated is 
proportional to (ey(v))*e lklyX ' 2 ~ Xl \ where A' is the quantum number of the particle that is removed 
by n z (x 2 ,t 2 ). As we shall see below, A' has to be equal to zero. This request is automatically 
guaranteed by another requirement, namely that the internal quantum numbers in the final state 
be exactly the same as those of the initial one. If this does not happen the final state is orthogonal 
to the initial one, and that particular configuration of quasiparticles does not contribute to the 
average (14). 

Let us now consider the case in which at time t = the point x\ — is located inbetween particle 
A p and A p+ i (sec figure (2)). Consider now a line connecting the points (0,0) and (0,ti) (this 
line does not correspond to any real particle trajectory). Define n' = N' + — N'_ as the number of 
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Figure 2. Propagation and scattering of quasiparticles in the semi-classical picture of the 0(3) 
non-linear sigma model. Lines drawn with different colours represent trajectories of different 
particles. The two states before and after the forward and backward time evolution must be 
identical, and this imposes constraints on the quantum numbers of the particles. Notice that in 
this example the particle with label p + 2 is removed at time 1 2 • 



intersections between this line and the other lines which correspond to particle trajectories, where 
N' + are intersections from the right and N'_ those from the left. By definition n' is given by: 

N 

n' = ^{eip-a^tiH-ep) 

V 

M 

= ^2{Q[Q-xi(h)] + Q[Q 

where the labels *- 1,2 ^ correspond to the trajectories of two quasiparticle originated at the same 
point, and N is the number of pairs N = M/2. O's are Heaviside step functions. Then at time 
t = t\ the operator h z (0, t\) creates a quasiparticle (L v = 1) with L z v — A„ = at x\ = with some 
momentum k. The adjoint particle will be created inbetween particles X' p+n , and A' +n / +1 . Suppose 
that n' > 0. The generalization of the following results to the case of negative n' is straightforward. 
The scattering 5-matrix takes on the exchange form (16), therefore during the motion particles 



x v }} 



xl{t 1 )]-2Q[Q-x u ]} 



(17) 
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only exchange their velocities while conserving their internal quantum numbers. This also means 
that the order of quantum number in each configuration from the left to the right never changes 
in time. The action of the operator n z (0,fi) on the initial sequence of quantum numbers may be 
written as: 

n z (0,ti) : {. . . , \p+ n ', Ap+ n '+i, . . .} ,^gN 
i-> {. . . , A p _|_„', Ao, A p+n '_|_i, . . .}. 

At time t = t 2 the operator n z {x,t 2 ) will remove the particle with label p + n 1 + n from the set 
of excitations, where n = N + — N_ is the difference between the number of intersections from the 
right and those from the left of the line connecting points (0,ii) and (x,t 2 ). This quantity is given 
by: 

N 
M 

= E i Q [ x ~ xl M] + @ l x - x lfa)] - 6[0 - xl(h)} - 0[O - a£(ti)]} . (19) 

This number tells us which particle is moving along this line at time t — t 2 , particle that will be 
removed by the action of h z (x,t 2 ): 

h z (x,t2) : {. . . , Ap+ n ', Ao, Ap+ n >+i, . . . , A p + n '+ n _i, Ap+ n '+ n , A p + n /+ n +i, . . .} ^q-j 

h-> {. . . , Ap+„', Ao, Ap+ n /+i, . . . , A p + n '+ n _i, Ap+ n '+ n +i, . . .}. 

Comparing the very last sequence of quantum numbers with the initial one (that at time t = 0), 
we end up with the following constraint on the set of quantum numbers: 



Ao = A p+ „' + i = A p+ „' + 2 • • • = Ap+„'+„ 



(21) 



In practice we have a sequence of n quantum numbers that must all be equal to A = 0, which 
starts at position p + n'. To identify the contribution of a particular configuration of quasiparticles 
we must consider the phase factors. Quasiparticles in the initial state generate a phase factor 
exp (— it 2 J2u e { v u)) under the action of exp(—it 2 H) (e(v) being the particle energy). This phase 
factor, however, completely cancels under the action of exp(it 2 H), except for the quasiparticlc 
added by the operator n z (0,ii). This gives a factor exp (— i(t 2 — t\)e(v)). Moreover every collision 
results in a sign change of the many-body wave function, but all these signs cancel under the forward 
and backward propagations, except those that are associated with collisions with the extra particle. 
These give an extra sign (—1) N++N - , which can be conveniently re-expressed as: 

Collecting all these pieces together we obtain the following expression for the general correlation 
function (7): 

C nlCTm 0ri, h; x 2 , t 2 ) = ^ (e (kye (k))c-^~ x i)-«W**-ti)) d kj x R( Xl , tl ;x 2 ,t 2 ), (23) 

where eo(fc) is the probability creation amplitude and R(x\, t\; x 2 , t 2 ) is the relaxation function, 
which is given by: 

+oo +oo 

R(x 1 ,t 1 ;x 2 ,t 2 ) = (J2 E (-^n'^l-A,!:^-] (24) 

n— — oo n'= — oo 

^ ,A p+n , +1 5A p+n , +1 ,A p+n , +2 • • • <5 A p+ „, + „_ 1 ,A p+ „, + „ ){x„,k„,\„}, ( 2 5) 
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/ 



where the sums with indices v\ and u 2 are given by equations (17) and (19) respectively. The (— 1)™ 
appearing in the previous equation is exactly that introduced in (22) . From the previous average we 
can immediately see that only those configurations which have a sequence of n quantum numbers all 
equal to A = contribute to the correlator (7). We now face the problem of computing an average 
of the form (0)i Xvi k vt \ u \, where O is a generic local operator, and the distribution probability takes 
into account all possible initial conditions. Solving this problem will enable us to compute the 
correlator (14). Due to the assumption of equation (15) we can write the average as: 

L/2 ^.-.^E E ••• E n^,A,)d({*^,M)- (26) 

aW.bW a( 2 ),f,( 2 > a (M), b (M)V=l 

fei>0 fe 2 >0 k M >0 

We will consider separately the cases in which (n',n) are even and/or odd numbers. 

• (n' = even and n = even): when n' is an even integer the adjoint particle is created inbetween 
two different quasiparticle pairs, thus we have a sequence of n/2 quasiparticle pairs which are 
all forced to have their quantum numbers equal to A = 0. The average over the set of quantum 
numbers can be written as: 

M 

= E E ••■ E n^(M 

a(i),b(i> a< 2 ),6< 2 ) a (M) b (M) „=1 [<t < ) 

X [*0.\ + „' + i^ p+ »'+i'Vn'+2 ' - - ^p+»'+»-i'W'+»_ ' 

where we have introduced the notation P a ' b (k v ) = P(k v ,X l/ ) (see figure (3)). 
It is worth noticing that the labels in the Kronecker deltas in the previous equation refer to the 
quasiparticles, while the labels in the summations and in the integrals refer to quasiparticle 
pairs. In order not to mix up these two different notations we recognize that in this case 
p + n' is an even number, therefore particle with label p + n' + 1 belongs to the 2: y L + 1 pair, 
and so on. From this pair to the right we have a sequence of n/2 quasiparticle pairs all with 
quantum numbers A = 0. In addition the probability distribution factorizes into single-pair 
probabilities. Let us remember that the average is taken over all possible initial conditions, 
therefore k v specifies the value of the momentum of the particle at time t = + e, where e is a 
small positive quantity. Roughly speaking k v is the value of the momentum of a particle right 
after time t = and before the first scattering process. Equation (27) can be easily evaluated 
to give: 

PittK}) = P(fci)...p(fc £ w)p 00 (fc £ ^ +1 )...P 00 (^ +5 ) 

' ' ( 2 8) 

xf(^ +f+1 )...P(M, 

where P(k) = J2a,b pab (k)- 

• (n' = even and n = odd): also in this case the adjoint particle is created inbetween two different 
pairs, but when n is odd the result of equation (27) becomes: 

p 2 «M) = P{k 1 )...p{k^)p m (k E ^ + ^...p w (k P _ ± ^ + ^) 

s „ 1 ✓ 

(29) 



E P° X {k e ^ + ^ +1 )p(k e ^ + ^ +2 )...P(k M ). 



A=0,±1 



Semi-classical theory for quantum quenches in the 0(3) non-linear sigma model 



11 




Figure 3. Computation of Pi({fc„}): an example with n = 6. Red trajectories represent those 
particles whose quantum number must be equal to Aq = 0. 



In practice the last particle of the sequence of A„ = belongs to a pair whose right-side partner 
can carry all possible quantum numbers. 

• (n' = odd and n = even): let us consider the case in which n' is an odd integer. In this case the 
adjoint particle is created within a quasiparticle pair, and the average over quantum numbers 
gives (n ^ 0): 



P 3 ({K}) = P(k 1 )...p(k P+n ,- l ) pX ° ( fc I 1 

2 A=0,±1 2 



X P QQ 



(30) 



J-J- — Itimcs 

^ \ ^ .-,1 1 \ / , 



E p0X ( fc ^- +t+ i) p ( fc ^ +t + ») • • • p ^- 



A=0,±1 



The case n = simply gives a sequence of M different P(fc„)'s. When n' is odd and n even 
in the sequence of quantum numbers with A„ = the first and the last one belong to pairs 
whose partner has no constraint on its quantum number A. This is the reason for the factors 
J2\ P ox (k) and J2\ P xo (k) in the previous expression. 
• (n' = odd and n — odd): in this case the average over the set of quantum numbers gives: 

P 4 ({M) - P(kl)...P(k p+ n>-l ) PX " ( fc * + " 2 '-l | 1 

2 A=0,±1 



X P QQ 



(31) 



•^4 — times 

x p(k ^-i , »-r |2 ) ...P(k M ). 
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Adjoint 
particle 




Figure 4. Computation of P2({k v }): an example with n ■ 



We are now ready to write down the global expression for the quantum number average: 

M 



ftp 



■ E 

a (M) b (M) v= l 
1 + (-1)"' 



(M<W p+ „, + /A p+ „, + 1 ,A p+ „, +2 



! + (-!)" 
2 

l + (-l) f 



1 - (-1)" 
2 

1 -(-!)" 



P2({K}) 



(32) 



This expression is a function of (n, n') and the set of quantum momenta of a particular quasiparticlc 
configuration. It is worth noticing that n and n' enter the expression for P$, i = 1, . . . , 4 in different 
ways. While n tells us how long the sequence is, n' gives us information about the starting point 
of the sequence and can easily be absorbed into the definition of p. As we shall see, for the actual 
computation of the relaxation function R the starting position of the sequence is not important, 
what matters is only the length of the sequence itself. We must now compute the average of 
the previous quantity over all possible initial positions and momenta. Namely we must take the 
following average: 



L / 2 dx\ dxM 



-L/2 



L 



+ 00 

E 

fci>0 k M >0 n= — oo 



E - - E E (-!)"*».&, [•••] E ^ 



(33) 



where the quantity in the squared brackets is given by equation (32). We can split this average into 
four pieces, let us call them A, B, C and D, one for each Pi({k u }) respectively, where i = 1, . . . , 4. 
We shall do the computation of the term which contains Pi({k v }) in some detail; the other cases 
are straightforward modifications. 

Let us compute explicity A, the first addend of the right side of equation (32), where the first 
summation over n' gets immediately canceled out by the corresponding Kronecker delta because nf 
itself enters the function Pi({fc y }) trivially, therefore we have: 

+oo 

E-E E c-in,Ej-] 

fci>0 k M >0 n=—oo 



A = 



L l' 1 dx\ dxu 



L 



L / 2 fei>0 fcj 



(34) 



Pi{{K}) 
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where ^2 [■ ■ ■] is specified by (17). A can itself be splitted into two pieces, say Ax and A 2 , which 
correspond to the two contributions to the integration of the terms into the bracket, namely: 

1 f-D^i 

2 + — 2 >A 1+ A 2 . (35) 

Let us start by computing the contribution of A\. By making use of the following integral 
representation of the Kronecker delta 

Ki:A.-] = ^r dte^-^A-]), (36) 



we can write A as: 



/ nn e in* l + (-!)" y y f L/2 dx x dx M c -^£,j...] ^ 
7t 27T — ' 2 — — J_t./o L L '' 

71 — — OO 



fei>0 fe M >o L / 2 

Now we want to compute the spatial integral, remembering that the sum 53r/ 2 [- ■ •] i s specified in 
equation (19). We begin by considering x > 0, t% > tj and X < v max (t 2 — ti), where u maa; is the 
maximal velocity of the excitations of the model (in the appendix B we will show the general result 
for arbitrary times and distances). We will refer to this situation as the within-the-light-cone case. 
The spatial integral factorizes in a straightforward way, thus we can write: 

f-L/2 7 j M 

r dxu e -^ E „ 2 [...] = JJ {j _ _ cosmVkJ x] 

J — L 1 i 



(38) 



-V2 - - „=1 

-^OKt - ar] - x (isin0 + 2sin 2 (0/2)) 9[z - v K t]Q[v k T - x] 

- {^f^{l + cos(20) -isin(20) - cos0 + isin0) + | (isin(2^») + 2sin 2 0)) Q[x - v kv T]\ , 

where Vk v is the particle velocity defined by = , t = t-2—t\ and T = t 2 + t\. Let us notice 
the following definition of the probability P(fc): 

p (*) = ^ M = X;/*, (39) 

fc>0 
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where f k is the occupation number of the fc-mode. Furthermore, the following identities can be 
useful in taking the thermodynamic limit: 

fc>0 _ k>0 JQ 27T k>0 _ fc>0 JQ 27T . . 

/• - /.» M ' /• ~ ./'/. ^ 

fc>0 fc>0 

Let us now compute the average over the momenta of the previous result, taking the limit L, M — > oo 
(M/L = p fixed): 

M 

fci>0 k M >0 i/=l 

where Pi({fci,}) is given by (28) and the {. . .} are the right-hand terms of equation (38). This 
expression can be simply evaluated to yield: 

M 

lim ]T Pi({M)I[(-) 

M.L-too ' — ' - L - L 
{fc„>0} f=l 

f f 00 dk f°° dk 

= exp <^ -t(l - cos0) / — /feVfceh^-x] -2xsin 2 (0/2) / — f k &[x - v k t]Q[v k T - x] 

I Jo * Jo 7r 

x exp ^ -T(l + cos(2^>) - isin(2</>) - 2cos0 + 2zsin0) / — / fe 6[x - v fc T] 
I Jo 2?r 

\ f°° dk 

x exp < — zxsin<^ / — (0[i>fct — x] + 0[x — v k t]Q[v k T — x]) 



(42) 



7T 



x exp^ -x(2sin 2 </>-Hsin(20)) / — f k Q[x - v k T]\(P m )\ n \' 2 , 

Jo 2tt J 

where P 00 = J* °° |^P 00 (fc) and P 00 (fc) is now a probability density. We can now take the limit 
when T is very large (x/T <C v max ), assuming that T J °° ^f k v k 8[x — v k T] — > 0, so the previous 
expression simplifies, giving: 



f [°° dk f°° dk 

exp <^ - cos0) / — fkV k @[v k t - x] - 2xsin 2 (0/2) / — / fe 6[x-i>fct] 
I Jo 7T Jo 7r 



x exp<j -ixsm^l —f k \(P 00 )^^ 



(43) 
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Before computing the sum over n it is convenient to obtain the analogy of equation (43) for the 
case A 2 , and see which contributions survive in the limit of T very large . For A 2 we start with: 



(44) 



fei>0 k M >0 L / 2 

then by repeating the same steps we end up with: 

M 

JTooE- E^«m)II(-) 

fci>0 fe M >0 ^=1 

(f°° dk f°° dk \ (ak\ 

-(T -tcoscj))j — f k v k - ix sin <j>j —f k <d[v k t-x]\ \ b > 

xcxp (-itsmjj™ ^.f kVk e[x-v k t}^j (p00)l«l/2. 

We immediately see that this contribution is exponentially decreasing as T becomes large, therefore 
we neglect it in this limit. We have already computed all the ingredients we need to do the sum 
over n, recalling that only even values of n give non-zero contribution, as can be seen from (37). 
Performing the summation, the global contribution from term A to the relaxation function (24) is 
given by: 

A = L l ™J £ (-d-^m ( 1+( - 2 1)E " 'V fi(W) 

e -t(l-cos0)/ o °° ^f k v k S[v k t-x] c -2xsin 2 (0/2)/~ &f k e[x-v k t] (46) 



00 dfc 



1 - (P 00 ) 2 



x cos (zsinH°° f /*) F _ 2p0 o cos(2 ^ + (P oo )2 > 

which in the large T limit depends only on t. Repeating all the steps we did for the term containing 
Pi({fc„}), we can then compute the contributions from the terms containing P2({fc;,}), Pz{{k u }) 
and P^({k v }) (namely P, C, and D). Here we list the results of these computations (always in the 
limit x/T < v max ): 

b = T ^ ( £ (-irs^.., ( 1+ (- 2 ^ i-(-D> 2({M) 

' {x„,k u } 

[ e -t(l-co S <p)f~ ff k v k e[v k t-x] c -2xsin 2 (<f>/2)f °°!ff k elx-v k t] (47) 



L,M-s-oo 

/+7T 



and 

+00 



\n— — 00 



/. 



{x„,fc„} 



"^g-tfl-cos^/^f/k^eKH c -2 2 ;sm 2 (0/2)/ o 00 e[x-«*t] (48) 
47T 



x cos fx sin 4 r * M 1 + 2c°s(2fl[(£ A P°*) 2 - P°°] - 2P°°(E A P°*) 2 + (P 00 ) 2 
xcos(xsm0j o ^/ fe j 1-2POOcos(20) + (POO)2 
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where we have assumed P 0A = P A0 , and finally 



\ n — — on ^ 



i-(-i) r 

2 2 



\n— — oo 



+7T 



e 



(l-cos«/ "f/ t » t 6K(-i] e -2xsin 2 (4>/2)/ °° e[s-i> fc t] (49) 

47T " 

/ ■ . r°°«ifc^ 2COS(/.(P 00 -1) ^ n0 A 

xcos(xsm0 /o fA) 1 _ 2 pooJ (2 ^ ) + ( poo )2 i: fOA - 

Plugging all the pieces together, A+B+C+D, we end up with the expression (11) for the relaxation 
function : 

/IT / roo \ 

— exp y-t(i -cos0)y — f k v k e[v k t- x]J 

/ f°° dk \ ( [°° dk \ 

x exp ( — x(l — cos(f>) J — fk Q[x — Vkt]j cos ( xsincf) J — fk 1 (50) 

1 + 2 cos(<ft) (P°° -1)Q + cos(20) (Q 2 - P°°) - P°°Q 2 
1 - 2P 00 cos(20) + (P 00 ) 2 ' 

where Q = X)a=±io^ , ° A - ^ would be useful to have an obvious comparison between this semi- 
classical result and a direct GGE computation of the same correlator for the case of the 0(3) 
non-linear sigma model, as was done for the transverse field Ising model and the XY chain in 
transverse field. The new approach to these kinds of problems introduced in [57] , [58] and [59] may 
be a missing tool in this respect. 

Formula (50) can be tested in several ways, for instance by choosing a thermally-populated initial 
state, with fk oc e~ l3k , and the same probability for each quantum number to appear. This check 
can be done both analytically and numerically. In the former case, starting from expression (50) 
and plugging in a thermal distribution for fk and the quantum numbers, we find the same universal 
analytical result of Reference [56] , whereas in the latter case we have performed a numerical average 
in the same spirit as that in Reference [60]. The results are shown in figure (7), which indicates a 
perfect agreement between the analitycal and numerical prediction. 



5. The dynamics of the transverse field Ising model 

5.1. The paramagnetic phase 

Another important check of formula (50) is represented by the transverse field Ising chain limit 
(TFIC). Briefly, we now focus on the dynamics after a quantum quench in the TFIC, whose 
Hamiltonian reads: 

h = -j ( E«+i - &X>n . ( 51 ) 

\i=l i=\ / 

where erf are the Pauli matrices at site i, J > is the energy scale, h the transverse magnetic field, 
and we impose free boundary conditions (for the moment we assume the system to be large but 
finite). The goal is to determine the dynamical order-parameter two-point function: 

C l ^{x,T,t) = ^\al 2 {t 2 )al i {t 1 m, (52) 
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Figure 7. The relaxation function R(x,t) for = Ae~@ k (where A = 0.1 and /3 = 0.8 have 
been used in this case) and p a > b = 1/9 for any choice of quantum numbers. Space and time are 
measured in unit of § and r respectively. Data points are numerical results and the solid line is 
the analytical prediction for the two-point function of a thermal initial state, see Eq. (50). 



where x = X2 — %x, T = and t = ti — t\. The action of the operator o~f(t) on the paramagnetic 

ground state is akin to that of the operator h z in the 0(3) non-linear sigma model, that is, it either 
creates an excitation or it destroys one already present. For this reason the semi-classical approach 
developed in the previous section is easily generalizable. For this model the K matrix in the def- 

where the 6's are the Bogoliubov 



■ tan 



inition of \ip) is known and equal to K{k) — 

angles defined in (55) before and after the quench respectively. In order to obtain (52), without 
rederiving everything from the beginning, one should formally replace P a ' b = l/(q — l) 2 and then 
take the limit q — > 2. Otherwise we can start again from (32), and take into account that for 
the TFIC excitations do not have internal quantum numbers (this is true both in the ordered and 
disordered phases). This feature greatly simplifies the algebra of the derivation, and we end up with 
the following general expression for the correlator of equation (52), valid in the thermodynamic limit: 



C Isin z(x,T,t) = C^(x,t) exp 
x exp 



rw fife 

-2t I — fkVk®[vkt-x 



IT I —f k v k Q[x - v k T] I exp 



r dk \ 
-2xj^ —f k e[x-Vkt}0[vkT-x}j, 



Semi-classical theory for quantum quenches in the 0(3) non-linear sigma model 



18 



Infinite Size Extrapolation ■ lt =5, h-2, !-30, T-16/3 




Infinite Sir,e Extrapolation - h -5, h-2, 1-30, T-J6/3 




Figure 8. The real and imaginary parts of the non-equal-time two point function after a quantum 
quench in the disordered phase, from hcj = 5 to h = 2. By definition t = t/tf and T = T/tt. 
Data points represent the numerical data extrapolated in the thermodynamic limit (L — > oo), 
while the solid line is given by equation (53). If we look at the inset in the plot of the real part, 
we notice that the semi-classical formula does not capture well the behaviour of the correlator 
outside of the light-cone (the same plot for the imaginary part shows a better agreement, even if 
at very small values of t numerical errors do not allow us to compare analytical result with the 
numerical ones). 



where C^™q (x,t) is given by: 

Crp_Q^X^ ^) OC 



dk Q~i e Wt+ikx 



2tt 



e(k) 



(54) 



where again as long as we work on the lattice the energy-momentum relation is e(k) = [(h— cos(fc)) 2 + 
sin 2 (/c)] 1 / 2 . Once you take the proper scaling limit a — > 0, then e(k) becomes relativistic and you can 
replace the general prefactor with a Bessel function Kq . Comparisons between theory and numerical 
data are shown in figures (8) and (9). We note that in expression (53) T is arbitrary (with the 
constraint T > t > 0), and not necessarily much longer than the Fermi time tf — x/v max , as was 
in (50). For larger quenches the semiclassical prediction, with ff. equal to \K(k)\ 2 /{\ + \K(k)\ 2 ) is 
not accurate. Indeed one has to take into account that for arbitrary values of the magnetic field 
h > 1 generic excitations are no longer single spin-flips (even if this is not the only source of error, 
see below) . In general their shape is given by a superposition of states with an arbitrary number 
of spin-flips, with coefficients depending on the value of h itself. Only when h 3> 1 arc single 
excitations well-approximated by single-spin flips, as you treat them in the semi-classical approach. 
A first correction to the semiclassical results corresponds to substituting — > —1/2 log | cos(Afc)|, 
where is the Bogoliubov angle given by: 

hoh — (ho + h) cos(fc) + 1 



cos(A fe ) = 



(55) 



e ho (k)e h (k) 

This substitution, which follows from asymptotically exact techniques (see [40]), increases the 
agreement between the theoretical predictions and numerical data, expecially close to the quantum 
critical point h = 1 where excitations are no longer localized objects. However in some plots 
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Figure 9. The real and imaginary parts of the non-equal-time two point function after a quantum 
quench in the disordered phase, from ho = 3 to h = 1.2. By definition t = t/tr and T/tj 3> 1. 
Data points represent the numerical data extrapolated in the thermodynamic limit (L — > oo), 
while the solid line is given by equation (53). The same considerations made for (8) hold true 
here. In addition we notice that as this quench is closer to the critical point, the agreement is 
only qualitative at small t, despite the substitution /j. — > — 1/2 log | cos(Ak)|. To improve the 
agreement we should correct the expression for the pre- factor of equation (53), following the exact 
results of reference [61] . 



mismatch between theory and numerics is still present, and we believe that it is due to the 
hypothesis of particle- number conservation of the pure semi-classical approach (i.e., the number of 
quasiparticles is conserved during the dynamical evolution of any configuration). If we considered 
the possibility of creating or destroying particles during the time evolution (which is realized by 
an operator insertion), we would get time-dependent corrections to the prefactor (54). Roughly 
speaking this would mean going beyond the leading order in the semi-classical approximation, a 
topic that is not addressed in detail in this paper. In the recent article [61] the dynamic correlators 
after a quantum quench were studied by using the form-factor technique, finding that the prefactor 
of expression (53) can be written as: 

dk c ikx T Ai 

t +2 i tan(-^)cos( efe (T))sgn(x- e ' fc t) . 



2 



(56) 



Outside the light-cone (when x > v max t) the first contribution in (56) is exponentially small, whereas 
the second one behaves as a power law. This observation explains in turn why the pure semi-classical 
approximation typically fails when applied to this regime. In particular the semi-classics is not able 
to capture the behaviour of the equal-time correlator, that represents the extreme out-of-the light- 
cone case. 



5.2. The ferromagnetic phase 

Formula (50) cannot be directly used to derive the correlator of the Ising model in the ferromagnetic 
phase. In this case excitations are no longer spin-flips, in fact when h < 1 these arc domain walls 



Semi-classical theory for quantum quenches in the 0(3) non-linear sigma model 



20 




Figure 10. Left panel (inset): real part of the non-equal-time (equal-time) two point function 
in the ordered phase. Again the numerical data are extrapolated in the thermodynamic limit 
(L — > oo). In this phase, unlike in the paramagnetic one, the agreement between the semi-classics 
and the numerics works well also for the equal-time correlator. Right panel: spatial correlation 
at equal times. 



(kinks). Following the general approach of F. Igloi and H. Reiger [37, 39] we can easily obtain 
the two-point correlation function for arbitrary value of times and distances. The key ideas of 
this method have already been extensively discussed in their works (even if these authors did not 
consider explicitly the case of dynamical correlators), here we just show the final result: 

C Ising (x, T, t) = C(h , h) R(x, T, t) (57) 

where the relaxation function in the ferromegnetic phase is equal to that in the paramagnetic phases 
while the multiplicative constant, as introduced in [42, 43] reads as: 

C(h ,h)= 1 - hh0 + ^ 1 -^ 1 -^. (58) 

In addition it is worth noticing that the semi-classical two-point correlator in the ordered phase of 
the Ising model is always a real function, as was the dominant contribution of the form-factor result 
at large (x,T) in [61]. In the ferromagnetic phase the equal-time correlators can also be described 
by the semi-classic approach, in contrast to the paramagnetic phase, as is shown in the inset of 
figure (10) 



6. Conclusions 



We have developed a semi-classical theory for the out-of-equilibrium quantum relaxation of the 0(3) 
non-linear sigma model, after having prepared the system in a coherent superposition of Cooper 
pairs, a structure that is in agreement with the integrability of the theory in the bulk. For such 
quenches we analyzed the two-point function of the order parameter fi z and argued that, in the long 
time limit its expression is given by formula (50), while for arbitrary times by equation (B.l). The 
method employed here is a generalization of that used for studying the finite-temperature behaviour 
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of a series of one-dimensional chains [49, 50, 56] . As was already observed in other integrable models, 
the long-time behaviour of this two-point function after a quantum quench is not thermal. For equal 
times (t 2 =ti^> x/v max ) the relaxation (50) reaches a stationary state, while for non-equal times 
(*2 ^ti) h is expressed as a function of the time difference t. These features are strong indications 
that the long-time limit of the two-point function (7) can be described by a statistical ensemble, 
and we expect it to be the GGE, being the model integrable. It would be interesting to have an 
independent and explicit result in the GGE framework, in order to check the general belief also 
for the 0(3) non- linear sigma model. The main difficulty here consists in finding an analytical 
expression for the lagrangian parameters of the GGE Hamiltonian, in order to explicitly carry out 
the statistical average. 

The semi-classical approach has also been applied to predict the dynamics of the order-parameter 
two-point function of the transverse field Ising chain, in both phases. In this case results for 
the auto-correlation and for the equal-time correlation are already present in literature [30, 37], 
the novelty here is the generalization of the semi-classical method to the case of different-times 
correlators. Yet, the exact results recently obtained in [61], have allowed us to explore the limit 
of the pure semi-classics. This technique works very well within the ferromagnetic phase (also for 
pretty large quenches), whilst it is in the paramagnetic phase where it shows its real limits: it fails 
in predicting the equal-time two-point correlation functions, and in general the agreement in the 
out-of-the-light-cone case (see the insets of figures (8) and (9)) is only qualitative. Furthermore in 
the proximity of the critical point, the method - which in general is supposed to work well only for 
small quenches - needs to be improved with the substitution f k — » —1/2 log | cos(A fe )|, which takes 
into account the real shape of the excitations. 

Nevertheless, our semi-classical theory can be applied straightforwardly to several other models, 
integrable and non-integrable, including the q-Potts model and the sine-Gordon theory, in order to 
compute correlators. These models will generate different combinatorial problems (the nature of 
the excitations is, of course, model-dependent), but the main ideas remain the same. 
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Appendix A. Occupation numbers and mode probabilities 

In this appendix we compute the occupation numbers and the probability distributions for each 
mode, starting from the definition of the squeezed coherent state (3) and the operator algebra of 
the operators (Z, Z^). Let us recall the definition of the initial state \ip): 




(A.l) 
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The Z operators obey the Zamolodchikov-Faddeev (ZF) algebra, namely: 
Z a {k x )Z h (k 2 ) - S c a t{k 1 ,k 2 )Z d {k 2 )Z c {k 1 ) = 

Zl(k{]Z\{k 2 ) - ScJ(k u k 2 )Z\{k 2 )Z\{k x ) = (A.2) 

Z a (kl)Zl(k2)-S^(k 1 ,k2)Zl(k 2 )Z c {k 1 ) = SabSih - k 2 ). 

Intuitively this means that the exchange of two quasiparticle is realized by the two-particle scattering 
matrix S(ki,k 2 ). In the present case we assume the form (16) for the S matrix, and thus the ZF 
algebra becomes: 

Z a {k x )Z h {k 2 ) + Zai^Zbih) = 

Z\{k{]Z\{k 2 ) + Z\(k 2 )Z\{k{) = (A.3) 

z a (h)zl(k 2 ) = [-Ec^MZM + sih -k 2 )]s ab . 

The ground state |0) satisfies Z a (k)\0) = 0, Va Vfc. From equations (A.2) it is clear that this algebra 
encodes the fact that the order of quantum number is preserved in time. The main quantities we 
want to compute are the occupation numbers f^ b , which are defined by: 

(^\zl(-k)zl(k)z b (k)z a (-kM b 

m = fk ■ ( 4) 

This quantity is exactly the number of quasiparticle pairs of kind (a, b) with momenta (— k, k) (see 
figure (3)), and it is directly connected to the probabilities P xx we defined in the body of the paper. 
The idea is to treat the if-matrix as an expansion parameter, but working from the beginning in 
the thermodynamic limit the squeezed state does not have a good expansion. Divergences appear 
in the expansion terms, in the form of squared Dirac delta-functions. In general there are two 
different methods of regularizing these divergences: one directly regulates the integral expressions 
in the infinite volume whereas the other operates through subtracting divergences in a large, finite 
volume. These two techniques were proposed and compared in [63]. In the follow we shall use the 
second method. 

Let us start by analysing the expansion of the denominator of (A. 4). First Taylor expand the 
squezeed state as: 

|V) = (1 + Y J K ab {k)Zl(-k)zl(k) + i(£ K«\k)Z\{-k)zl{k)f + . . .)|0>, (A.5) 

fc>0 fc>0 

a.b a, b 

therefore we have: 

(V#) = l + Z 1 +Z 2 + ..., (A.6) 

where 

fc>0 £>0 k>Q 

a, b c,d a,b 

As expected this term is proportional to the volume L. the higher orders will be proportional to 
L 2 , L 3 et cetera. Let us consider the expansion of the numerator of (A. 4): 

^\zi(-0zl(0Zp(0z a (-^M) = w 1 + \w 2 + ..., (A.8) 

where the first term Wiis given by: 

^ = E E<°i^(*2)^(-*2)^(-o4(o^(o^a(-o^(-*i)^(*i)io> = i^coi 2 - (a. 9 ) 

fcl>0 fc 2 >° 

a,b c,d 
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This quantity is finite, also in the infinite volume limit L — > oo, whilst W 2 is proportional to L. We 
now need the expression of W% to see how these divergences cancel against the normalization of the 
boundary state. By definition we have: 

w * - E E E E<°i z «( fci )^(- fci )^( fc2 )^(- fc2 ) z «(-o4(o 

fc!>0 fc 2 > D fc 3>° fc 4> D (A 10) 

a, 6 c.c£ e.f g.h ' ' 

x^(0^(-0^(-fc3)4(A;3)^(-fc4)4(^)|0) 
After a bit of algebra we end up with the following result for W 2 : 

W 2 = 4|^(0| 2 Yl \K a \k)\ 2 + 2\K^(0\ 2 J2\ Kab ^)\ 2 , (A.11) 

fe>o a, 6 

a, 6 

where we immediately see that the first term is of order L. In particular the first term of (A. 10) 
can be written as AW\Z\, thus the first correction in equation (A. 4) is given by: 

m) = i + Zl + ■ ■ ■ = Wl + o{Wl) - (A - 12) 

We see that the divergent term indeed cancels against the normalization factor, and this happens 
order by order. We conclude that as long as K ah (k) is small we can calculate the matrix elements 
like (A. 4) by expanding the squeezed state. The final result can be written as: 



ff = \K ab (k)\ 2 +o(\K ah (k)\ 2 ) (A.13) 



which explains the name amplitudes for the matrix elements K ab (k). The momentum occupation 
numbers fk are given by: 

/* = E/£ 6 < (A-14) 

a, b 

while for the probability densities P ab (k) -which arc defined in the thermodynamic limit- we have: 



fk b _ 

dk 



pab ( k ) = /eg AU , (A.15) 



2tt 



•ab 

k 



a.b 

which makes the normalization condition J2 ab J^° ^P ab (k) = 1 explicit. The probability P ab is 
defined by P ab = J °° ^P ab (k). 



Appendix B. General expression of the relaxation function R(x,t,T) 

In this appendix we show the form of the relaxation function for the 0(3) non-linear sigma model 
when the times T and t are arbitrary, with the only constraint T > t > 0. This expression can be 
written as: 



R(x;t,T) = R 1 {x;t,T) + R 2 {x;t,T), 



(B.l) 
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where the fisrt addend reads as: 



dd> ( f°° dk 

— exp - cos 0) I — f k v k Q[v k t - x] 

-7T 27T V J Q 7T 

{!f 00 dk \ 
-T(l + cos 20-2 cos 0) I —fkVk @[x — v k T] 
Jo 2lT ' J 

f°° dk \ 
—x(l — cos 0) / — fk &[x — Vkt]Q[vkT — x] I 
Jo 7T ' / 
r°° dk \ 
-z(l - cos 20) / — / fc 6[a; - v k T] j 
f°° dk ° [°° dk 

xs\n2(j)J — fk®[x - v k T] + xsin<j) j —f k Q[v k T-x\ 



x cos 



r°° dk 1 
sin 20)/ —f k v k Q[x- v k T] \ 

rpoo _ un i r .r^(oA>\(r>2 _ pdo\ 



+ 2T(sin - sin \ 

1 + 2cos(0)(P°° - j)Q + cos(20)(Q 2 - P tf0 ) - P m Q 2 
1-2P 00 cos(20) + (P 00 ) 2 



(B.2) 



and the second term is given by 

R 2 (x;t,T) = J ^cxp^-2(T-icos0)^ ^f k v k B[v k T - x] 

( f°° dk \ 

x exp I -x(l - cos 20) / -fkQ[x-v k T]\ 

) 1 r°° dk \ 

x exp f [-T(l + cos 20) + 2t cos 0] / 7^fkV k Q[x - v k T] J 

f f°° dk ° f 00 dfc 

xcosj2xsin0/ — f k <d[v k t - x] + 2t sin <j) — f k v k @[x - v k t]Q[v k T - x] 

r°° dk r°° dk 1 

+ (-Tsin20 + tsin0) / —f k v k e[x-v k T]+xsm2(j) — f k G[x - v k T] \ 

P°°Q 2 - (P 00 ) 2 - cos O (20)( 7 Q 2 - P°°) 
X 1 - 2P 00 cos(20) + (P 00 ) 2 

This formula represents the leading order of the semi-classical approach to the n z — n z correlator in 
the 0(3) non-linear sigma model. From this expression all subcases can be derived straightforwardly. 



(B.3) 
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